Theoretical analysis of continuously driven dissipative solid-state qubits 
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We study a realistic model for driven qubits using the numerical solution of the Bloch-Redfield 
equation as well as analytical approximations using a high-frequency scheme. Unlike in idealized 
rotating- wave models suitable for NMR or quantum optics, we study a driving term which neither 
is orthogonal to the static term nor leaves the adiabatic energy value constant. We investigate the 
underlying dynamics and analyze the spectroscopy peaks obtained in recent experiments. We show, 
that unlike in the rotating- wave case, this system exhibits nonlinear driving effects. We study the 
width of spectroscopy peaks and show, how a full analysis of the parameters of the system can be 
performed by comparing the first and second resonance. We outline the limitations of the NMR 
linewidth formula at low temperature and show, that spectrocopic peaks experience a strong shift 
which goes much beyond the Bloch-Siegert shift of the Eigenfrequency. 



Coherent manipulation of quantum states is a well es- 
tablished technique in atomic and molecular physics. In 
these fields, one works with "clean" generic quantum sys- 
tems which can be very well decoupled from their envi- 
ronments. Moreover, it is possible to apply external fields 
in a way such that strong symmetry relations between the 
static and the time-dependent part of the Hamiltonian 
apply and the resulting dynamics is very simple and can 
be treated analytically. In solid-state systems, the situa- 
tion is different. Not only do they contain a macroscopic 
number of degrees of freedom which form a heat bath 
decohering the quantum states to be controlled, but also 
is the choice of controllable parameters much more re- 
stricted. A quantum-mechanical two state system (TSS) 
realized in a mesoscopic circuit can be identified with a 
(pseudo)spin, however, in that case the different com- 
ponents of the spin may correspond to physically dis- 
tinct observables such as e.g. magnetic flux and electric 
charge 1 . This naturally limits the possibilities of con- 
trolling arbitrary parameters of the pseudospin. Hence, 
in order to describe the direct control of quantum states 
in mesoscopic devices, concepts from NMR or quantum 
optics cannot be directly applied but have to be carefully 
adapted. In particular, as decoherence is usually rather 
strong in condensed matter systems, one can attempt to 
drive the system rather strongly in order to have the op- 
eration time for a quantum gate, usually set by the Rabi 
frequency, as short as possible. 

We concentrate on the case of a persistent cur- 
rent quantum bit 2-4 driven through the magnetic flux 
through the loop and damped predominantly by flux 
noise 5 with Gaussian statistics. This setup is accurately 
described by the driven 6 ' 7 spin-boson model 8 
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where e(t) = eo + scosQt and the oscillator bath is 

assumed to be ohmic with a spectral density J(ui) = 

2 

f Y,i ^b"^ - ^) = aue- u ^ c . The connection of J(u) 
to the setup parameters is detailed in 5 . The sta tic en- 
ergy splitting of the peudospin is v = Veg + A 2 . This 
model is also applicable to other Josephson qubits and 
other realizations 8 ' 9 . In particular, the strong driving 
regime we are going to elaborate on has recently been 
realized in several setups 10-12 . We study the effective 
dynamics of the pseudospin having traced out the bath 
in the limit of weak damping, a <C 1, which is appro- 
priate for quantum computation. This is done using the 
Bloch-Redfield equation 13 . The resulting equation is of 
Markovian form in the sense that it only contains the den- 
sity matrix at a single time, however, it is derived in such 
a way, that the free coherent evolution during the inter- 
action with the bath is fully taken into account such that 
the resulting equation is numerically equivalent to a fully 
non-Markovian path-integral scheme 7 ' 9 and only mem- 
ory terms beyond the Born approximation are dropped. 
The explicit form of the equations for this situation as 
well as the formulas for the rates correspond to those 
given in 7 . We compare our numerical results to analyti- 
cal formulae derived in the framework of a high-frequency 
approximation 14 which involves averaging over the driv- 
ing field and has nonetheless shown to give a good esti- 
mate for the system dynamics even close to resonances 7 . 

Initial experiments on quantum bits such as Ref. 3 do 
not monitor the real-time dynamics of the system as in 
Ref. 4, because the read-out is much slower than the de- 
coherence, i.e. the dephasing time is too short. In 
order to optimize the experimental setup, it is important 
to measure both and the relaxation time tr, even 
and in particular if they are insufficient. In the stan- 
dard NMR-case, this is done by studying the width of the 
resonance 15 . We will detail, a somewhat modified analy- 
sis can be performed for solid-state qubits and what are 
its limitations. We discuss both situations 3 ' 4 . Our re- 
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suits thus help to analyze the decoherence as observed in 
Refs. 3, 4, and outline the possibilities and limitations of 
driving the system in the nonlinear regime. We have nu- 
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FIG. 1. (a z )(t) at fixed frequency Q/2tt = 6.6GHz for dif- 
ferent bias points, a) Off-resonance, v/2-k — 2.9GHz, inco- 
herent decay towards thermal equilibrium; c) on-resonance, 
v /2ix — 6.6GHz, Rabi oscillations decaying on the scale of a 
dephasing time towards a dynamical equilibrium state; b) 
close to the resonance, v /2n — 6.4 GHz, combination of deco- 
herence and relaxation; d) short-time dynamics highlighting 
the fast oscillating component, see text. 

merically solved the driven Bloch-Redfield equation. The 
real-time dynamics is illustrated in fig. 1. 

The dynamics shows distinct features on different time 
scales. As expected, there are clear Rabi oscillations on 
the scale of the effective driving strength (see below). In 
quantum computing applications, these would be used 
for the implementation of a Hadamard gate. On top 
of this, there are fast components: The dominating one 
oscillates with the driving frequency, which originates in 
the fact that the driving is not perpendicular to the static 
field. A weaker one, which oscillates at twice the driving 
frequency, comes from the counter-rotating term perpen- 
dicular to the static field. These oscillations can lead to 
errors of the Hadamard gate. On a longer time scale, the 
Rabi oscillations decay. The time scales will be discussed 
later on. In general, if one is not exactly on resonance, 
these oscillations are combined with nonoscillatory de- 
cay, see fig. 1 a) and b). At very long times, the system 
assumes a quasistationary value P^. 

Corresponding to the situation of a spectroscopy ex- 
periment, we now turn to the analysis of the quasi- 
stationary state which is established after a long time 
t ^> T^iTRiUJft 1 . We compare our full numerical solu- 
tions with analytical expressions we have obtained from 
the high-frequency approximation of Refs. 7, 14. As a re- 
sult of this approach, the TSS is mapped onto a coupled 



ensemble of TSSs corresponding to the original system 
emitting or absorbing n photons from the driving field 
during the tunneling. The energy bias of these individual 
systems is e n = eo — nfl and the tunnel matrix element 

A n = AJ n (s/n) (2) 

where the J n are Bessel functions. At low driving fields, 
we can approximate A n = ^(s/2Q) n as we would expect 
from the expansion of a perturbation series in the driv- 
ing strength. The A n can hence be viewed as n-photon 
Rabi frequencies. This implies, that the usual single- 
photon frequency gets replaced by Ai ~ sA/v, which can 
be interpreted as only the projection of the driving field 
onto the direction in pseudospin space orthogonal to the 
static Hamiltonian. In order to obtain the solid curves in 
fig. 2 the secular equations for the Eigenfrequencies have 
been solved, taking into account an appropriate number 
of terms 16 . The dynamical two state systems are charac- 
terized by individual dynamical dephasing rates r<^ n and 
a common relaxation rate T r 7 . On the n-th resonance, 
can be very low, much lower than off-resonance, as 
can be seen in figure 1, and largely exceed the intrin- 
sic dephasing time. This has been observed in Refs. 1, 
4. Figure 2 shows numerical and analytical results for 
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FIG. 2. (<j z ) in the long-time limit at fixed frequency 
Q/2tt = 6.6GHz as a function of the energy bias eo for differ- 
ent values of the driving strength s/Q = 0.034,0.43,1.7,2.4 
(a to d). From b to d, nonlinear resonances can be identified, 
d already shows negative values at small positive eo, which 
can be identified as the coherent destruction of zero-photon 
tunneling. Further peaks occur at even higher bias. 

Poo = \im t ^oo(o'z) at a fixed frequency Q/2tt = 6.6GHz 
as a function of the energy bias eo- This corresponds to 
a realistic experimental situation 3 . In fig. 2 a), taken at 
weak driving field, only the regular resonance correspond- 
ing to the transition between the two Eigenstates driven 
by absorbing a single photon can be seen. At somewhat 
stronger driving, fig. 2 b), this peak grows wider and a 
second resonance appears, corresponding to the simulta- 
neous absorption of two photons. At higher fields, fig. 2 



c), these peaks grow and start to dominate over the back- 
ground. They also turn asymmetric. This trend culmi- 
nates in the situation shown in fig. 2d). In that case, 
does not grow to positive values at small positive eo, but 
it gets negative and then directly approaches the first 
resonance. The reason for this behavior can be identi- 
fied within the high-frequency approximation: The low- 
est order-tunnel frequency Ao = AJo(s/Q) vanishes at 
this particular driving strength. Indeed, comparing figs. 
2 a) and d) one can see, that the step which is at eo = in 
case a) is shifted to eo — Q in case d). This phenomenon, 
the coherent desctruction of tunneling 6 relies on destruc- 
tive interference of the dressed state 17 formed by the TSS 
and a cloud of photons from the driving field. This in- 
terpretation is supported by the dynamics of (a z )(t). As 
seen in fig. 3, which shows the dynamics at the degener- 
acy point for different driving strengths, the zero-photon 
tunneling is slowed down and brought to a standstill. If 
that strong driving can be applied to solid-state qubits, 
it would provide an alternative for controlling A by a 
cw microwave field instead of an additional magnetic flux 
as proposed in Ref. 2. At very weak driving, the peak 
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FIG. 3. (a z ) at strong driving with high frequency, 
Q/2tt = 6.6 GHz (where A/2tt = 660 MHz) By increasing the 
driving strength, the tunneling is slowed down and brought 
to a standstill. 

position corresponds to the qubit eigenfrequency, Q = v. 
This is not reliably predicted by the high-frequency ap- 
proximation. At stronger driving, the peak gets shifted. 
Closer inspection as in figure 4 shows, that this shift goes 
much beyond the usual Bloch-Siegert shift 6 of the dy- 
namical Eigenfrequency, in fact, one can show that the 
position of the 'peak in steady state and the Eigenfre- 
quency do not conincide. The former is given by balanc- 
ing of rates and it can be shown that in lowest order gets 
shifted by 16 ^e pea k — s 2 /8Q 2 whereas the Bloch-Siegert 
shift for our case is 5cbs — A 2 s 2 /(16^ 3 ). As a more 
general conclusion, already at modest not-too-weak driv- 
ing, the resonance positions do not necessarily reflect the 
Eigenfrequencies of the system. In fig. 5, the height of 
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FIG. 4. Shift of the spectroscopy peak as a function of the 
driving strength for uj/2tv = 6.6GHz and A/2tt=660 MHz. 
We compare to the usual Bloch-Siegert shift formula (dashed) 
and the formula derived in the text (solid). 
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FIG. 5. Heights of the two lowest resonances as a function 
of the driving strength s at is/2tt = 6.6 GHz, A/2tt = 660 
MHz, Q/2tt = 6.6 GHz The solid and dashed lines are extrap- 
olated NMR formulae. 



the two lowest order peaks is shown. It can be seen, that, 
from the low-driving side, they saturate as soon as their 
effective Rabi frequency A n exceeds 1 / ^T r T$. At very 
high driving, the peaks show an inversion of population. 

For the optimization of qubit setups on the way to 
coherent dynamics, it is important to characterize its co- 
herence properties from the spectroscopic data. In NMR, 
this is done from the linewidth given by 



sn 



(3) 



where uor is the Rabi frequency at resonance, which co- 
incides with the strength of the driving field. A gen- 
eralization of this formula to our case has to take into 
account low temperatures and the different driving situ- 
ation. Moreover, uor is usually not directly known to suf- 
ficient precision, because the driving strength depends on 



the attenuation of the applied fields on their way to the 
sample and the efficiency of the coupling 3 ' 5 . Our analysis 
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FIG. 6. Widths of the two lowest resonances as a function 
of the driving strength s at A/2tt = 660 MHz Q/2tv = 6.6 
GHz. The solid line corresponds to the extrapolated NMR 
formula discussed in the text. 



suggests the generalization of eq. 3 is given by 



SQ n = 2Jr t 
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(4) 



where 5fl n is the width (in frequency) of the n-photon 
resonance and A n is the effective Rabi frequency defined 
above. At low powers s < they are given by the rates 
from the undriven ohmic case 
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This result is confirmed by our numerical simulations fig. 
6. We can essentially identify three regimes: A saturation 
broadening regime at low powers, where 5Q n 2r^ 1 , a 

saturated regime, SQ n ~ 2A n y^r^/r^ and a nonlinear 
regime, where the numerical curve deviates from eq. 4 
due to the fact, that the high Rabi frequency shifts the 
relevant energy scales and modifies the timescales given 
in eq. 5. Note, that in this regime, the general curve of 
Poo is greatly deformed (see fig. 2) and the width of a 
peak becomes ambiguous. 

This result allows to measure essentially all interesting 
parameters of the system experimentally. By extrapo- 
lating the level separation at the degeneracy point (as 
it was done in 3), one obtains A. By tracking the reso- 
nance positions at weak driving, one can evaluate eo as 
a function of the external control parameter (in 3 this 
would be the magnetic flux). By driving in the saturated 
regime, the widths of the first and second peak become, 
according to eq. 4 £Qi/ 2 = 2A 1 / 2 \AW r 0> hence by tak- 
ing their ratio we find the effective driving strength from 



A2/A1 = Ji(s/Q)/ Jo(s/Q) ~ s/2 and by tracking the 
slope of the first resonance we find the ratio tv/t^. Fi- 
nally, examining the saturation broadening regime of the 
first resonance gives the absolute value of r^. 

In conclusion, we have numerically and analytically 
analyzed the spin-boson system, which e.g. represents a 
SQUID-qubit, in the weak damping regime, driven by 
continuous fields. As compared to the more familiar sit- 
uation in NMR, this system is both different in the char- 
acter of the driving and the low temperature governing 
the dissipation. We have shown, that the key features 
of this system, Rabi-oscillations, and saturation of the 
linewidth, persist qualitatively as has been experimen- 
tally confirmed 4 . They are however altered on a quanti- 
tative level, such as an unanticipatedly strong shift of the 
position of the resonance peak, and also supplemented 
by new phenomena such as higher-harmonics generation, 
oscillations of (a z ) on the scale of the driving field, and 
coherent desctruction of tunneling. We have finally out- 
lined a scheme how to determine all relevant parame- 
ters (tunnel splitting, energy dispersion, driving strength, 
dephasing and relaxation time) of a quantum bit solely 
through spectroscopy. 
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